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Abstract. The aim of this paper is to examine the biogeographical distribution of the littoral and sub-littoral opistho- 
branch gastropods from Brazil and nearby areas. On the basis of published literature and personal data, a zoogeographic 
study was undertaken and the Brazilian region was compared with the Caribbean and Argentinean regions. The Brazili¬ 
an littoral has been divided into six zones, based in the oceanographic features established by Castro & Miranda (1998). 
466 species belonging to the orders Cephalaspidea, Anaspidea, Sacoglossa, Notaspidea and Nudibranchia from the Ca¬ 
ribbean to the Argentinean region were considered in this study. The number of opisthobranch species is highest in Ca¬ 
ribbean areas. Along the Brazilian littoral. South Brazilian Bight is the region with highest richness, while Orinoco-Ama- 
pa and Amazon shelf are the areas with lower species numbers. The similarity analysis shows that some possible geo¬ 
graphic barriers act to the distribution of the opisthobranchs. For each locality considered in the study, the percentage of 
species that extends northward is higher than southward. 
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1. INTRODUCTION 

Faunistic inventories give fundamental information for 
many basic and applied scientific disciplines, such as ecol¬ 
ogy and biogeography (Stork & Samways 1995). An es¬ 
sential tool to establish inventories is taxonomy, which 
supplies a reference system for the biodiversity (Bisby 
1995). The level of faunistic and taxonomical knowledge 
from different regions varies considerably, which raises 
a problem for the accomplishment of biogeographic stud¬ 
ies. Nevertheless, the elaboration of this type of studies 
offers insights into possible models of distribution of a tax¬ 
onomic group throughout a more or less extensive area. 

Opisthobranch gastropods are well represented in most 
marine habitats from equatorial to polar regions. Ernst and 
Eveline Marcus studied this group for more than thirty 
years in West Atlantic temperate and warmer waters. Re¬ 
cently, other authors have provided new contributions to 
the knowledge of the Brazilian opisthobranchs (Troncoso 
et al. 1998; Garcia et al. 2002; Garcia & Troncoso 2003, 
2004; Padula & AbsalAo 2005; Pola et al. 2005; 
Dominguez et al. 2006 a, b; Valdes et al. 2006). The ac¬ 
tual level of knowledge of the opisthobranchs from 


Brazilian and Caribbean regions (i.e. Marcus 1977; Mar¬ 
cus & Marcus 1967a; Thompson 1977, 1980; Espinosa 
& Ortea 2001; Rosenberg 2005; Valdes et al. 2006) per¬ 
mits some biogeographical considerations; however, fu¬ 
ture studies along some areas of Brazilian coast are nec¬ 
essary to obtain a more complete knowledge of this fau¬ 
na. The work of Marcus and Marcus was focused around 
Sao Paulo and Rio de Janeiro areas. The necessity of more 
faunistic studies in Brazil is obvious when comparing the 
number of species cited along Brazilian coasts (280 
species, after Marcus 1977) with, for example, the Iber¬ 
ian Peninsula (523 species, after Cervera et al. 2004). 

Here, we present a study on the diversity of the Opistho¬ 
branchs along the Brazilian coasts, using as biogeograph¬ 
ic areas the six zones defined by Castro & Miranda 
(1998). These Brazilian zones are also compared with 
Caribbean and Argentinean regions (i.e. Muniain 1997; 
Schrodl 1999; Rosenberg 2005). 
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2. MATERIAL AND METHODS 


2.1. Biogeographica! areas 

This research was conducted by comparing littoral and 
sub-littoral Opisthobranch fauna from the Brazilian 
shores with those from the Caribbean and Argentinean bio- 
geographie regions. Species checklists were compiled by 
combining data from bibliographical sources and person¬ 
al observations (included in the references). The total num¬ 
ber of species considered for this study is 466, belonging 
to the orders Cephalaspidea, Anaspidea, Saeoglossa, No- 
taspidea and Nudibranehia (Table 1). 

In this paper, we adopted the six oeeanographie zones de¬ 
fined in Castro & Miranda (1998) for Brazil (Fig. 1). 
The geographical limits and features of these areas are 
shown in Table 2. 



Fig. 1. Limits of the areas considered along Brazilian coasts, 
based on Castro & Miranda (1998) and nearby zones based 
on Briggs (1974). The number of species considered for each 
region is indicated. 


Within the Caribbean region we have considered the 
provinces defined in Briggs (1974) (Fig. 1): 

CC-CR, extends from Cape Canaveral to Cabo Romano 

CR-CRo, extends from Cape Romano to Cape Rojo 

CRo-O, extends from Cape Rojo to the mouth of the 
Orinoco River 

O-A, extends from Orinoco River to Amapa 
ABC, Includes Aruba, Bonaire and Curaeao 
WI West Indies (WI) 

Bermudas Is (BER). 

In addition, we consider the Argentinean region (ARG) 
between the borders of Brazil to Uruguay up to 43-44°S 
(Chubut). 


2.2. Community analysis and measurement of biodi¬ 
versity 

In order to examine diversity within the opisthobranch 
communities, data by region were subjected to a multi¬ 
variate analysis using the Bray Curtis similarity measure 
and non-metric Multidimensional Scaling Ordination 
(MDS). 

The Bray-Curtis index (Bray & Curtis 1957) was cho¬ 
sen because it does not consider double absences in its cal¬ 
culations. The results were then graphically described 
using dendrograms with the UPGMA (unweighted pair- 
group methods using arithmetic averages) aggregation al¬ 
gorithm (Sneath & Sokal 1973). The ordination analy¬ 
ses were carried out by means of an MDS (non metric mul¬ 
tidimensional sealing program) based on the similarity ma¬ 
trix between stations. 

For two different station groups a requirement is to iden¬ 
tify which species account for the observed assemblage 
difference (Clarke & Gorley, 2001). The SIMPER rou¬ 
tine was used to identify taxa that greatly contributed to 
differentiate station groups. The software used was 
P.R.I.M.E.R. (Plymouth Routines in Multivariate Ecolo¬ 
gical Research) version 5.2.8. for Windows. 

3. RESULTS 

466 species of opisthobranehs were found to oecur in the 
Western Atlantic Ocean from Cape Canaveral to the Ar¬ 
gentinean province (the Atlantic Magellanic region is not 
included). Table 3 shows the number of speeies by order 
or suborder for each zoogeographic area. The number of 
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opisthobranch species for each area varies remarkably. WI 
(West Indies) has the highest number of species (250 
species), although, in general, the number is also high in 
the Caribbean areas Cro-0 and CC-CR. Along the 
Brazilian zones, the number of species is lower and there 
is also a notable difference between the different zones. 
SBB is the zone with the highest number of species while 
AS is the zone with the lowest number (122 and 12 
species, respectively). The number of species from the Ar¬ 
gentinean region (ARG) is moderately low. 


3.1. Faunal affinities 

The cluster analysis using the six Brazilian zones shows 
a first division in which AS separate from the remaining 
zones (Fig. 2A). A second division separates SBS, and fi¬ 
nally. the rest of the areas splits into two groups; a group 
including the north-eastern Brazilian zones EBS and NBS, 
and the other including ACR and SBB. A two-dimension¬ 
al representation of the analysis MDS shows this same 
grouping pattern (Fig. 2B). 




Fig. 2. Cluster classification (A) and MDS ordination (B) of 
the biogeographic areas based on the presence-absence of spe¬ 
cies in the areas from Brazil. 
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Fig. 2. Cluster classification (A) and MDS ordination (B) of 
the biogeographic areas based on the presence-absence of spe¬ 
cies in all areas included in this study. Black triangles, Brazili¬ 
an regions; white triangles. Caribbean provinces; black circle, 
Argentinean region. 


Including all geographical areas from Cape Canaveral and 
Bennudas to Chubut (Argentinean region), the cluster 
analysis and MDS show a first division in which 
Orinoco-Amapa (O-A), the Brazilian areas SBS and AS 
and the Argentinean region (ARG) separate from the re¬ 
maining regions (Fig. 3). Those separate in two groups. 
One includes the north-eastern Brazilian zones EBS and 
NBS, while the other includes the remaining zones. The 
latter divides itself into two subgroups, one including the 
Brazilian regions SBB and ACR, and the other subgroup 
including all Caribbean provinces. 

The SIMPER analysis identified four distinct groups. In 
group 1 (SBS, O-A, ARG, AS), the following taxa con¬ 
tribute to the similarity (up to 85%); Cephalaspidea 
(Acteocina bicientata, Volvulellcipersimilis, Acteocina can- 
clei , Acteon pelecais) and Notaspidea (Pleurobraiichaea 
inconspicna). Group 2 (WI, CC-CR, CRo-O, CR-CRo, 
ABC. Ber) is mainly characterised by species of the or- 
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ders Cephalaspidea ( Volvulellci, Bulla, Haminoea, Hvdati- 
na, Micromelo ), Notaspidea ( Umbraculum, Pleuro- 
branchus ), Anaspidea ( Stylocheilus, Bursatellci, Aplysia , 
Bosellia ), Sacoglossa ( Cylimlrobulla, Oxynoe, Tridachia, 
Elysici) Dendronotacca (Scyllaea) and Nudibranchia 
(Spurilla, Chromodoris). Group 3 (SBB, ACR) is charac¬ 
terised by the presence of species of the order Nudi¬ 
branchia, mainly Doridina (Dendrodoris krebsii, Cadlinci 
rumia, Tyrinna evelinae, Diaulula greeleyi, Discodoris 
evelincie, Chromodoris clenchi ), Arminina (Annina) and 
Acolidina (Spurilla neapolitana, Phidianci lynceus, Fla- 
bellina engeli, Glaucus atlanticus). Group 4 (EBS. NBS) 
is determined (up to a cumulative 90%) by the presence 
of the orders Nudibranchia (Doto divae, Diaulula gree¬ 
leyi ), Notaspidea (Pleurobranchaea inconspicua, Umbrac¬ 
ulum umbraculum, Berthellina quadridens ), Anaspidea 
(Phyllaplysia engeli. Sty 'loch ei Ins stri at us, Aplysia dacty- 
lomela, Aplysia parvula, Bursatella leachii ), Saccoglos- 
sa (Cylimlrobulla beauii, Oxynoe antillarum, Elvsia tuca ), 
and Cephalaspidea (Philine sagra, Hydatina vesicaria, Mi¬ 
cromelo undatus, Haminoea elegans, Clielidonura petra) 

In terms of dissimilarity, the species Spurilla neapolitana, 
Scyllaea pelagica, Pleurobranchus areolatus, Tridachia 
crispata, Aplysia dactylomela, Aplysia fasciata, Oxynoe 
antillarum contribute greatly to differentiate groups 1 and 
2. Group 1 differed from group 3 due to Phidiana lynceus, 
Spurilla neapolitana, Flabellina engeli, Facelina coenda, 
Anteaeolidiella indie a, Berghia benteva. Taringa telopia, 
Tyrinna evelinae, Siraius ilo, Jorunna spazzola, Okenia 
zoobotiyon, 11 all ax a apefae, Doris verrucosa, Discodoris 
evelinae, Chromodoris neona, Dendrodoris krebsii, Diaul¬ 


ula greeleyi, Cadlina rumia, Chromodoris clenchi, 
Bertliella agassizii, Berthella stellata, Aplysia dacty'- 
lomela, Aplysia fasciata, Oxynoe antillarum, Ascobulla ul- 
la, Navanax aenigmaticus. Group 1 differed from group 
4 due to Doto divae, Diaulula greeleyi, Phyllaplysia en¬ 
geli, Stylocheilus striatus, Berthellina quadridens, Aplysia 
dactylomela, Aplysia parvula, Oxynoe antillarum, Ely si a 
tuca, Micromelo undatus, Chelidonura petra, Atys ri- 
iseanus. Differences between group 2 and group 4 were 
mainly due to Tridachia crispata, Aplysia fasciata, Che¬ 
lidonura petra. The group 2 differed from group 3 due to 
Facelina coenda, Berghia benteva, Scyllaea pelagica, Sir¬ 
aius ilo, Hallaxa apefae, Chromodoris neona, Tridachia 
crispata, Philine mera, Cylindrobulla beauii, Acteon pele- 
cciis. And finally the species Phidiana lynceus, Flabelli¬ 
na engeli, Glaucus atlanticus, Facelina coenda, Annina 
tnuelleri, Anteaeolidiella indica, Berghia benteva, Doto di¬ 
vae, Taringa telopia, Tyrinna evelinae, Siraius ilo, Oke¬ 
nia zoobotiyon, Hallaxa apefae, Chromodoris neona, 
Dendrodoris krebsii, Cadlina rumia explained most of the 
dissimilarity between groups 3 and 4. 

Figure 4 shows the number and percentage of endemic 
species for each biogeographic area. Those species con¬ 
sidered as endemic have been cited only at one zoogeo¬ 
graphic region from the area of study. Geographic distri¬ 
bution along other biogeographic regions was not consid¬ 
ered for this study. The level of endemism varies notably 
along the different zoogeographic areas. The highest val¬ 
ue was found in the Argentinean region (68.6%). This high 
value is due to a southward distribution to the Magellan 
region of the fauna from this area. The Brazilian zones 
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Fig. 4. Number and percentage of endemic species for each faunistic area. 
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have percentages between 0 % in Amazon Shelf (AS) and 
East Brazilian Shelf (EBS), and 28.7 % in South Brazil¬ 
ian Bight (SBB). Among the Caribbean provinces, the per¬ 
centage of endemism varies between 7.1 %, at Aruba, 
Bonaire and Curacao (ABC), and 22 % in the West Indi¬ 
an region. 

Table 4 shows the percentage of species from each region 
present in other localities. The percentage is increasing, 
when extended northward, and decreasing, when extend¬ 
ed southward. 


4. DISCUSSION 

An aspect to be considered in this study is that the data 
on species distribution and the level of knowledge of the 
communities vary across the geographical regions. Thus, 
the results presented here must be considered as tentative. 

The number of opisthobranch species between biogeo¬ 
graphic areas varies remarkably being higher in the 
Caribbean: West Indies (Wl: 250 speeies). Cape Rojo- 
Orinoco (CRo-O: 178) and Cape Canaveral-Cape Romano 
(CC-CR: 173). The south Brazilian Bight (SBB) is the 
zone in Brazil with the highest richness (122 species), fol¬ 
lowed by the north Brazilian Shelf (NBS: 75) and the 
Abrolhos-Campos Region (ACR: 62). The areas with the 
lowest species richness are Oricono-Amapa (O-A) and 
Amazon Shelf (AS), both with 12 species, followed by 
South Brazilian Shelf (SBS) with 18 species. This could 
be related to oceanographic conditions, such as the effect 
of the Amazon River in the AS zone, and with a differ¬ 
ence in sampling effort; the south Brazilian Bight is the 
zone where Ernst and Eveline Marcus did many of their 
collections. 

The first division observed in the similarity analysis in¬ 
cludes four geographical regions with features that seem 
to act as barriers to the distribution of opisthobranchs. The 
coasts from Orinoco-Amapa and Amazon Shelf are influ¬ 
enced by the Amazon and Orinoco rivers whose plumes 
spread north-westward for more than 1,000 km into the 
North Atlantic (Castro & Miranda 1998). These areas 
are characterised by soft bottom, turbid waters and fresh¬ 
water runoff and they have been recognised as barriers to 
the dispersal of corals (Cox & Moore 2000), rocky shore 
gastropods (Vermeij 1978) and shallow water reef fish¬ 
es (Gilbert 1972). The scarce opisthobranch fauna from 
Orinoco-Amapa and Amazon Shelf and the composition 
of species present in these areas are likely related to en¬ 
vironmental characteristics of the region. Species present 
in these areas are generally cephalaspideans and arrni- 
nacean, which frequently live in sand or mud. 


A similar situation is found on the South Brazilian Shelf, 
where the low species numbers may be related to the ef¬ 
fect of the Patos Lagoon River plume with an annual mean 
discharge of about 2000 mV 1 (Marques et al. 2006). This 
area is influenced during the winter by Subantarctic wa¬ 
ter. 

The coastal area of the Argentinean biogeographic 
province is a transition zone characterised by processes 
of mixing and instability of the water masses. I his 
province includes geographical features like the Rio de la 
Plata. The estuaries influence the primary and secondary 
production in the area and consequently, the distribution 
of species. Moreover, there is an interaction off the coast 
between the Malvinas current flowing on the slope from 
the south with cold Subantarctic waters rich in nutrients 
and the Brazilian current, with temperatures higher than 
20° C and salinity over 36.0 ppt. This determines the pres¬ 
ence of eurythermal and euryhaline species (Boschi 2000). 
The fauna of opisthobranchs is formed mainly by nudi- 
branchs, which have their northern distribution limit at the 
border between the Argentinean region and the South 
Brazilian Shelf, extending southwards to Subantarctic re¬ 
gions. 

Cluster analyses show two Brazilian groups, composed of 
NBS-EBS and ACR-SBB. In general, these groups coin¬ 
cide geographically with those indicated by Floeter et al. 
(2001) to the reef-fish fauna of the Brazilian coast. These 
authors considered several regions like the South and 
South-eastern coastal reefs, from the Guarapari islands to 
Santa Catarina (areas included in our analysis as the group 
ACR-SBB), and the North-eastern coast, extending from 
the Manuel Luis reefs to Abrolhos Archipelago (areas in¬ 
cluded in our analysis as the group NBS-EBS). 

The southern and south-eastern coastal reefs, cited by 
Floeter et al. (2001) show lower mean annual water tem¬ 
perature, relatively higher primary production and a large 
shelf width. The reef-fish fauna living in this area appears 
to be the richest of Brazil, due to the mix of tropical and 
subtropical elements. The area is subjected to a relative¬ 
ly intense seasonal upwelling promoted by the South At¬ 
lantic Central Water, bringing low-temperature (<18°C) 
and nutrient-rich waters close to the coastline (Ekau & 
Knoppers 1999). Floeter et al. (2001) stated that a con¬ 
siderable number of Caribbean reef fishes found in this 
region are absent from the north-eastern sites. We found 
a similar pattern in opisthobranchs. The cluster analysis 
groups ACR-SBB closer to the Caribbean provinces than 
to other Brazilian zones. On the other hand, ACR and SBB 
are the Brazilian areas with the highest richness in opistho¬ 
branchs. This could be related with environmental features 
(as is observed for reef-fishes; Floeter et al. 2001). Nev¬ 
ertheless, in addition to environmental factors, differences 



208 


Francisco J. Garcia et al.: Biogeography of Brazilian Opisthobranchia 


in the richness of opisthobranehs for each area may de¬ 
pend upon other factors, such as discrepancies in sampling 
effort. The fauna from Brazil is better known in the South 
Brazilian Bight and Abrolhos-Campos zone, where Eve¬ 
line and Ernst Marcus conducted research for over 30 
years. 

The north-eastern region as is described by Floeter et al. 
(2001) for fishes, which nearly overlap with the group 
NBS-EBS, is characterised by its relatively warm waters, 
a weak seasonal signal and a small vertical temperature 
gradient; the circulation is influenced northward by the 
North Brazilian Current, and southward by the Brazilian 
Current (Castro & Miranda 1998). The reef formation 
consists of eoralline algal erusts over a rocky substrate, 
hennatypie and fire corals, as well as sponges (Floeter 
et al. 2001). The narrow and open shelf is an oligotroph- 
ie system almost entirely covered by carbonate sediments 
due to little freshwater input and the coast is influenced 
by the South Equatorial Current (Knoppers et al. 1999). 
In these areas herbivorous Saeoglossa and Anaspidea are 
more abundant. This trend of an increase in abundance in 
herbivores towards the tropical zone was previously ob¬ 
served in the Brazilian reef fishes (Ferrfira et al. 2004). 

With regard to the distribution of opisthobraneh species, 
it can be noted that for each locality considered in this 
study, the percentage of species that extend northward is 
higher than southward. Several authors discussed the 
Southern Caribbean as a centre of origin of species (see 
Briggs 2006) from which the species have been penetrat¬ 
ing northward into Florida and Bermudas and southward 
into Brazilian waters (Rocha 2003). Concerning Opistho- 
branehia, more intensive faunistie studies along the South 
American Atlantic coast are needed, before this hypoth¬ 
esis can be applied to understand opisthobraneh diversi¬ 
ty in that region. 
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Table 2. Limits and features of the areas considered along Brazilian coasts, based on Castro & Miranda (1998). 
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Table 3. Number of species by order or suborder for each area considered. 






